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We consider a massless scalar field propagating in a weakly curved spacetime whose metric is a 
solution to the linearized Einstein field equations. The spacetime is assumed to be stationary and 
asymptotically flat, but no other symmetries are imposed — the spacetime can rotate and deviate 
strongly from spherical symmetry. We prove that the late-time behavior of the scalar field is identical 
to what it would be in a spherically-symmetric spacetime: it decays in time according to an inverse 
power-law, with a power determined by the angular profile of the initial wave packet (Price falloff 
theorem). The field's late-time dynamics is insensitive to the nonspherical aspects of the metric, 
and it is governed entirely by the spacetime's total gravitational mass; other multipole moments, 
and in particular the spacetime's total angular momentum, do not enter in the description of the 
field's late-time behavior. This extended formulation of Price's falloff theorem appears to be at odds 
with previous studies of radiative decay in the spacetime of a Kerr black hole. We show, however, 
that the contradiction is only apparent, and that it is largely an artifact of the Boyer-Lindquist 
coordinates adopted in these studies. 



I. INTRODUCTION AND SUMMARY 

This paper is concerned with the late-time behavior of 
radiative fields in stationary curved spacetimes, a topic 
that has been the subject of vigorous investigation in the 
last three decades. Here we take the spacetime to be that 
of an isolated body of mass M, and we take the metric to 
be stationary and asymptotically flat. We imagine that 
the radiative field (a classical, massless field of integer 
spin that propagates on the fixed background spacetime) 
is given suitable initial data on a hypersurface t — 0, 
and we examine its evolution at times t much larger than 
M. Under these conditions, it is found that the field's 
late-time behavior is characterized by an inverse power- 
law decay. This property was first discovered by Richard 
Price [1], and it has been extensively studied since. This 
paper makes an additional contribution to this topic's 
vast literature by generalizing the formulation of Price's 
falloff theorem: we show that it holds under less restric- 
tive conditions than are normally assumed. Although 
our considerations here will be restricted to a scalar field, 
much of what we shall say should apply equally well to an 
electromagnetic field, or to a (linearized) gravitational- 
perturbation field. 



Radiative falloff in spherical spacetimes 

Excluding the more recent work on radiative decay in 
the Kerr spacetime, which we review below, all previous 
studies of the late-time dynamics of massless fields were 
restricted to spherically-symmetric spacetimes. Here 
the field admits a decomposition in terms of spheri- 
cal harmonics, and the field equations reduce to a one- 



dimensional wave equation with an effective potential. 
For example, the evolution of a massless scalar field ^' in 
Schwarzschild spacetime is governed by the reduced wave 
equation 
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where r* = r + 2M\n{r/2M + 1) is the usual tor- 
toise coordinate, V{r) = (1 - 2M/r)[l{l + l)/r^ + 
2M/r^] is the effective potential, and ipi{t,r) is a 
dimensionally-reduced field from which the scalar field 
can be reconstructed. The relationship is r^(t, r, 9, cj)) — 

J2im"'i-rn'4"-(^y'^)^r'i^^^)^ where aim are amplitudes de- 
termined by the initial data, and Y™{6, </>) are the stan- 
dard spherical-harmonic functions. With a small change 
to the effective potential, 2M/r^ going to 2{1 — s'^)M/r^, 
the wave equation of Eq. (1.1) can be related to a massless 
field of (integer) spin s [1, 2], and this generalized poten- 
tial can be shown to lead to a universal late-time behav- 
ior for ^pi; conclusions based on the scalar-field equation 
therefore apply to fields of higher spin. 

In his original treatment [1] , Price dealt with the par- 
tial differential equation of Eq. (1.1) and solved it as a 
perturbative expansion in powers of M. This method 
was extended in Refs. [3, 4]. These works reveal that for 
times much larger than r (which can be taken to be much 
larger than M), ipi decays in time according to ^-(2^+3)^ 
assuming that the field's initial data has compact sup- 
port and is not time-symmetric. This limit corresponds 
to approaching * = 00 keeping r finite; this gives the 
field's behavior near future timelike infinity. Another in- 
teresting limit corresponds to setting v = t + r* =00 
keeping u = t — r* finite; this gives the field's behavior at 
future null infinity. In this limit ipi behaves as it^('+2) 
The behavior of the field on the event horizon was also 
determined [3, 4], but we shall not discuss this here. 

In an alternative strategy. Leaver [5] performed a 
Fourier transform (in t) of Eq. (1.1) and examined the re- 
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suiting ordinary differential equation (in r* ) . He showed 
that the inverse power-law decay of the field is associated 
with a branch-cut discontinuity of the solutions in the 
complex frequency plane, and he reproduced the correct 
late-time behavior with a low-frequency approximation 
{Mw <C 1, where w denotes frequency) to the differential 
equation. The low-frequency regime is naturally related 
to the field's late-time dynamics {t » M), and to the 
asymptotic behavior of the effective potential (r ^ M). 
This method was generalized in Refs. [6-8]. 

The work of Ching et al. [6] is especially notewor- 
thy, because it contributes a very simple heuristic un- 
derstanding of the inverse power-law decay. These au- 
thors show that the late-time field at the spatial posi- 
tion r is produced by the part of the initial wave packet 
that propagates outward to a distant point r" (such that 
r" is much larger than r), scatters off a small amount 
of spacetime curvature there, and makes its way back 
to r after a time t 2r". They show further that if 
5V{r) = V{r) - l{l + l)/r*2 is the part of the effec- 
tive potential that can be directly associated with space- 
time curvature, then the late-time field is proportional 
to 6V{t/2)-. this holds for arbitrary potentials V{r). The 
inverse power-law nature of the effective potential nat- 
urally gives rise to an inverse power-law decay for the 
field. 

The work reviewed above reveals that the late-time be- 
havior of a radiative field is governed only by the asymp- 
totic form of the metric at large distances from the grav- 
itating body. In this sense, the radiative falloff brings 
information about the spacetime's asymptotic structure 
to an observer located much closer to the gravitating 
body. On the other hand, the late-time dynamics of a 
massless field reveals nothing of the nature of the central 
object, which might be extremely diffuse or extremely 
compact. While the body's internal structure certainly 
affects the field at early times [7, 9-17], the radiative 
falloff depends only on the asymptotic conditions. The 
robustness of the inverse power-law decay is well demon- 
strated in Refs. [18-20], in which numerical simulations 
of the collapse of a self-gravitating scalar field were car- 
ried out; the decay is seen at late times whether or not 
the collapse produces a black hole. That the late-time 
dynamics of a massless field is directly coupled to the 
spacetime's asymptotic structure is well illustrated in 
Refs. [21-27], in which different asymptotic conditions 
are imposed, and corresponding changes in late-time be- 
havior observed. Changes also occur if the simple mass- 
less fields considered here are replaced by more exotic 
fields; this was explored in Refs. [28-32]. 



New formulation of Price's falloff theorem 

Our previous conclusion, that the late-time dynamics 
of a massless field is governed only by the spacetime's 

asymptotic structure, and is independent of the details 
of the metric near the center of mass, is based on calcu- 



lations that assume spherical symmetry. But this obser- 
vation suggests that the assiimption is not a severe re- 
striction, as any stationary, asymptotically-flat spacetime 
becomes approximately spherical at large distances from 
the gravitating body. We would therefore expect that the 
picture of radiative falloff painted previously should not 
change qualitatively if the assumption of spherical sym- 
metry were to be removed. Our main purpose with this 
paper is to prove that this expectation is warranted: We 
investigate the radiative falloff of a massless scalar field 
in a nonspherical spacetime, and show that the same pic- 
ture emerges. 

The observation also suggests that in order to study 
the late-time dynamics of a massless field, it should be 
sufficient to examine its propagation in a spacetime that 
is only weakly curved, since this condition is necessarily 
met in the asymptotic region. This, then, defines our ap- 
proach: We look at the late-time evolution of a massless 
scalar field ^ in a stationary, asymptotically-flat space- 
time whose metric satisfies the Einstein field equations 
linearized about flat spacetime. 

The metric is expressed as 

ds^ = -{l+2^)dt^+{l-2^){dx'^+dy^+dz^)+S{A-dx)dt, 

(1-2) 

in terms of a scalar potential <^{x) and a vector potential 
A{x), in which x = {x,y,z) denotes a quasi-Cartesian 
frame adapted to the underlying flat spacetime. The po- 
tentials do not possess any (spatial) symmetries, but it 
is assumed that they do not depend on time t. The mat- 
ter distribution can thus be arbitrarily distorted, but it 
is assumed to be stationary and bounded by a sphere of 
radius TZ. The spacetime is therefore asymptotically flat, 
and the potentials outside the matter distribution are 
characterized by a (potentially infinite) number of mass 
and current multipole moments. Details regarding the 
metric and the potentials are presented in Sec. II. De- 
tails regarding the equation satisfied by the scalar field 
are given in Sec. III. Throughout the paper we work con- 
sistently to first order in the gravitational potentials, and 
we use units such that G = c = 1. 

From the quasi-Cartesian coordinates {x, y, z) we form 
quasi-spherical coordinates (r, 9,(f)) by invoking the usual 
relations x = r sin 9 cos (j), y — r sin 9 sin 0, and z = 
rcosO. We then give the scalar field initial data 

^{o,x') = c{r')Yr{e'A'), 

(1.3) 

^{0,x') = C{r')Yr{6'A') 

on the hypersurface t = Q. (The symbol x' designates 

a point on the hypersurface.) Here, C(r') represents the 
radial profile of the initial wave packet, and C'(r') gives 
the profile of the field's velocity. [These functions are in- 
dependent and the overdot should not suggest that C{r') 
is the time derivative of C(r'); it is just a convenient 
mnemonic that reminds us that while C(r') is associ- 
ated with the field's initial configuration, C{r') is associ- 
ated with the field's initial velocity.] We assume that the 
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field's initial data has compact support: the fimctions 
C{r') and C{r') must vanish outside of a sphere of ra- 
dius C. This is an essential assumption: as we shall see, 
the methods developed in this paper cannot be directly 
applied to initial data that do not satisfy this property. 
As Eq. (1.3) indicates, we assume! that the field's initial 
angular profile is described by a spherical-harmonic func- 
tion of degree I and azimuthal index m. This choice is 
largely arbitrary, because the absence of spherical sym- 
metry in the background spacetime denies any special 
status to the spherical harmonics. They constitute, nev- 
ertheless, a convenient basis in which to decompose the 
field's initial angular profile. The additional choice of as- 
signing the same spherical harmonic to both ^'(0, x') and 
^'(0, x') is also arbitrary, and it is made for the sake of 
simplicity. 

The initial data of Eq. (1.3) is evolved forward in time 
in the weakly curved spacetime of Eq. (1.2). This is done 
by means of a retarded Green's function G{t,x;t' ,x') 
that we compute in Sec. IV. In that section we prove one 
of the central results of this paper: Despite the absence 
of spherical symmetry, for times t > t' + r + 3r' + 2TZ 
the Green's function is sensitive only to the spherically- 
symmetric aspects of the metric; it is entirely oblivious 
to the vector potential A and the nonspherical part of 
the scalar potential We show that for these times, the 
Green's function is given by 



G{t,x;t',x') 



8M{t - t') 



\x — X 
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(1.4) 



where |a; — a;'| is the Euclidean distance between x and 
x', and M is the total gravitational mass of the space- 
time. Equation (1.4) indicates very clearly that only 
the monopole moment of the scalar potential enters in 
the description of the late-time Green's hmction; higher 
multipole moments play absolutely no role. It should be 
noted that in Sec. IV, the Green's function is calculated 
perturbatively through first order in the gravitational po- 
tentials; Eq. (1.4) is therefore valid up to terms quadratic 
in $ and A. 

This remarkable property of the Green's function guar- 
antees that in the weakly curved spacetime of Eq. (1.2), 
the late-time behavior of the scalar field will be the same 
whether or not the spacetime is spherically symmetric. 
While the early-time behavior brings out essential differ- 
ences, wave propagation at late times is insensitive to any 
nonspherical aspect of the spacetime. The spherical pic- 
ture of radiative falloff painted previously is therefore not 
just a good approximation, it is virtually exact, so long 
as we do not probe the nonlinear aspects of the problem 
(see below for a clarification of this qualifying statement). 

By integrating the initial data of Eq. (1.3) against the 
Green's function of Eq. (1.4), we find in Sec. V that the 
scalar field's behavior near future timelike infinity {t r) 
is given by 



*(i»r) ~ 4M(-1) 



i-i-i 



(2/ + 2)!! 
(2^ + 1)!! 
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f2l+3 



(1.5) 



up to a fractional correction of order jt^ . The quan- 
tities Ci and Ci are moments of the field's initial radial 
profile; for example, C; = J C{r'){r'y^'^ dr' . In Sec. V 
we also obtain the field's behavior at future null infinity. 
To get this we introduce the retarded time coordinate 
u = t — r, the advanced time coordinate v = t + r, and 
we take the limit f — » oo. After multiplying the field by 
r to produce a nonzero result, we find 



r^/{v = oo) ~ 2M(-1) 



i+i 



(2/ + 1)!! 



(1.6) 



up to a fractional correction of order C? jv? . These re- 
sults agree precisely with those obtained previously for 
Schwarzschild spacetime [3, 5]. In Sec. V we show that 
both relations can be obtained by taking a different 
limit of a more general description according to which 
cx r'+it/(t2 - r2)'+2 oc (1 - 7i/w)'+i(l + ?i/?.)/«/+2. 
The transition from a ^"(^'+3) behavior near future time- 
like infinity to a u~('+2^ behavior at future null infinity 
is thus naturally explained. To the best of our knowl- 
edge, the interpolating formula of Eq. (5.15), below, has 
not appeared before in the literature; that Eqs. (1.5) and 
(1.6) can be obtained by taking two different limits of a 
more general description appears to be a new result. 

Our new statement of Price's falloff theorem is consid- 
erably more general than previous formulations, which 
were all restricted to spherically-symmetric spacetimes. 
We prove here that in a stationary, asymptotically-flat, 
and weakly-curved spacetime without spatial symme- 
tries, the late-time behavior of a scalar field is identical 
to what it would be in a spherically-symmetric space- 
time. We find that the late-time dynamics is insensitive 
to the nonspherical aspects of the metric, and that it is 
governed entirely by the spacetime's total gravitational 
mass; other multipole moments, and in particular the 
spacetime's total angular momentum, do not enter in the 
description of the field's late-time behavior. This conclu- 
sion is compatible with our previous intuition: Because 
the radiative falloff is governed by the asymptotic form 
of the gravitational field, and because the field must be- 
come increasingly spherical at large distances from the 
gravitating body, the spherical description should be re- 
liable. What we show here, in effect, is that the spherical 
description is virtually exact, and that it is even more 
reliable than could have been expected. 



Radiative falloff in Kerr 

This conclusion appears to be at odds with previous 

studies of radiative decay in the spacetime of a Kerr black 
hole. This topic has been the subject of both analytical 



4 



[33-38] and numerical [39-41] investigation. For the pur- 
pose of this discussion we rely mostly on the analytical 
description of the late-time behavior of a scalar field pro- 
vided by Hod [37] . The most reliable numerical work also 
focused on the radiative fallofF of a scalar field, and it was 
carried out by Krivan [41]. 

The late-time dynamics of a massless scalar field in 
the Kerr spacetime appears to be much more compli- 
cated than what is suggested by our preceding discus- 
sion. Supposing once more that the field's initial data is 
a pure spherical- harmonic mode V™, the field's behavior 
near future timelike infinity is found to be given by [37] 

r y™/t2i+3 l = morl = m + l, 
* oc I y^^^/t'+^+i l-m>2 (even), (1.7) 
[ y^JViA'^"^^ l-m>2 (odd), 

in apparent conflict with the result stated in Eq. (1.5). 
The Kerr results imply that an initial mode of degree I 
does not necessarily survive at late times. In general it 
is converted into a mode of degree m or m -t- 1; this is 
the mode of smallest degree compatible with the initial 
value of m (which docs not change because of the axisym- 
metry of the Kerr spacetime) and with the parity of the 
initial mode (which is conserved during the evolution). 
Ignoring the fact that Eq. (1.7) does not agree with our 
previous predictions, that evolution in the Kerr space- 
time should produce mode conversion is not necessarily 
surprising: the spacetime is not spherically symmetric, 
and spherical-harmonic modes should not be expected to 
evolve independently. It is therefore expected that an ini- 
tial Y;"* mode should be converted to a mode if Z — m 
is even, or to a Y^+i mode if Z — m is odd. What is sur- 
prising, however, is that a final y,™ mode does not neces- 
sarily decay in time according to i-(2'"+3)^ which is what 
it should do according to our spherical-symmetry expec- 
tation; instead, the decay is given by ^-C+^+i) when 
I > m + 2. Similarly, a final Y™_^i mode does not de- 
cay according to i-(2"»+5)^ i,^^ goes instead as ^-('+'"+2) 
when I > m + 2. 

The behavior displayed in Eq. (1.7) is loosely confirmed 
by Krivan's numerical work [41]. Krivan shows that an 
initial mode is converted to a Yq mode during the 
evolution, and finds a decay described approximately by 
^-5.5^ This is much closer to Hod's prediction of 
than to the naive (spherical-symmetry based) expecta- 
tion of t~^. Although the agreement is far from perfect, 
it clearly suggests that Eq. (1.7) is closer to the truth 
than Eq. (1.5), even after allowing for mode conversion. 

The field's radiative falloff at future null infinity is even 
more complicated. Here Hod finds [37] 

( y;™/u'+2 l = morl = m+l, 

oc <^ (y^. ■ ■ ■ ^"2)/"'' l-m>2 (even), 

I (iT+i, n;"+3, ■ • • Yr^^)/u^ l-m>2 (odd), 

(1-8) 

where the notation indicates that in general, the angu- 
lar description of the late-time field involves a number of 



spherical-harmonic modes. The mode-conversion mech- 
anism is clearly more complicated here, and the results 
of Eq. (1.8) disagree with Eq. (1.6); they also seem to be 
incompatible with our spherical-symmetry expectations. 

Our new formulation of Price's falloff theorem makes 
an unambiguous statement: The late-time behavior of 
a scalar field in a stationary, asymptotically-flat, and 
weakly curved spacetime is governed only by the spher- 
ical aspects of the metric, and is not affected by any 
deviation from spherical symmetry. Provided that the 
strong-field aspects of the spacetime can truly be ignored, 
as we have argued, the theorem predicts that a scalar 
field in the Kerr spacetime should behave as predicted 
by Eq. (1.5) and (1.6). Instead we find the behavior of 
Eq. (1.7) and (1.8). Is there a contradiction between the 
two sets of results? 

The answer is no there is no contradiction. The 
key point is that the late-time dynamics of a scalar field 
in Kerr was examined in the Boyer-Lindquist coordi- 
nates, which give one of the simplest descriptions of the 
Kerr spacetime. But in these coordinates, the asymp- 
totic form of the Kerr metric does not match the met- 
ric of Eq. (1.2), whose relevant spherical part is given 
by ds2 ^ _ 2M/r) df + (I + 2M/r){dr^ + dO^ + 

siv? 6 d(f>^). Indeed, even in the limit M — > the Kerr 
metric does not coincide with the Minkowski metric ex- 
pressed in spherical coordinates (r,0,<j}). To produce a 
match it is necessary to introduce spheroidal coordinates 
{f,9,(p) defined by the relations x = ^/f'^ -\- sin ^ cos 0, 
y = \/ + sin 6 sin (/), and z = rcosO, in which the 
constant e measures the ellipticity of the coordinate sys- 
tem. Rewriting the Minkowski metric in terms of these 
coordinates produces agreement with the M — » limit of 
the Kerr metric, provided that e be identified with Kerr's 
rotation parameter: e = J/M, where J is the spacetime's 
total angular momentum. 

It is therefore the spheroidal coordinates, and not 
the original spherical coordinates, that constitute the 
appropriate weak-field limit of the Boyer-Lindquist co- 
ordinates. This implies that a meaningful comparison 
between the results of Eqs. (1.5), (1.6) and those of 
Eqs. (1.7), (1.8) must involve a transformation from the 
spherical coordinates (r, 6, cf)) — used to express the first 
set of results — to the spheroidal coordinates (f, 6, (j>) — 
implicitly used to express the second set of results. Stated 
differently, the initial data used in the context of Kerr has 
the form of C {f)Y™ {9 , (p) , and this is substantially dif- 
ferent from what is written in Eq. (1.3). Equations (1.7) 
and (1.8) refer to this different choice of initial conditions, 
and they do not agree with Eqs. (1.5) and (1.6) simply be- 
cause the evolutions proceed from different initial data. 
The two sets of results therefore represent very differ- 
ent situations, and it is meaningless to compare them 
directly. 

We elaborate this argument in the last two sections 
of the paper. In Sec. VI we retrace the computations 
of Sec. V, starting with initial conditions of the form of 
Eq. (1.3) but stated in terms of spheroidal coordinates. 
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We evolve this initial data with the Green's function of 
Eq. (1.4), also expressed in terms of spheroidal coordi- 
nates. Our results for the late-time behavior of the scalar 
field are in perfect agreement with Hod's results [37] dis- 
played in Eqs. (1.7) and (1.8). In Sec. VII we give a direct 
comparison between the descriptions of radiative falloff 
based on spheroidal and spherical coordinates, and show 
how they can be reconciled. The analysis presented in 
these two sections implies that the results of Eqs. (1.7) 
and (1.8) are largely an artifact of the Boyer-Lindquist 
coordinates, which obscure the effective spherical sym- 
metry of the problem. It also shows that there is no 
contradiction between Hod's results and our generalized 
statement of Price's falloff theorem: the higher multi- 
pole moments of the Kerr spacetime, including its an- 
gular momentum, play no role in the description of the 
field's late-time dynamics. 



Open questions 

Because our calculations agree completely with Hod's 
[37], they shed no light on the apparent discrepancy be- 
tween the analytical description of Eq. (1.7) and Krivan's 
numerical results [41]. (Recall that Krivan finds an ini- 
tial mode to decay approximately as t~^'^, while Hod 
predicts with an error of order i~^.) A possible ex- 
planation for this discrepancy might be that Krivan's nu- 
merical results are sensitive to subdominant terms in the 
expansion of \1/ in inverse powers of t; this would explain 
the slightly larger power. It is also possible that Kri- 
van's results are revealing a nonlinear effect of the sort 
discussed by Hod in Ref. [42] . More work will be needed 
to clarify this issue. 

It will also be important to better understand the 
limitations of our weak-field assumption. Consider the 
early-time evolution of a scalar field from t = to some 
time t = t\. During this epoch the field is sensitive to 
the nonspherical aspects of the metric, and the initial 
data of Eq. (1.3) gets converted into a superposition of 
(potentially many) spherical-harmonic modes. Of these, 
the most relevant for the field's late-time behavior is the 
mode of lowest degree in I, because it is this mode that 
will dominate at late times. (For example, an initial 
field should be expected to produce a Y^ mode at t = ti , 
and this should give rise to a field that decays as t~^, 
not t~^^ , at late times.) Isolating this mode, the field's 
new initial data at t = ii can once more be written as in 
Eq. (1.3), but with a potentially lower value oil, and with 
radial functions C(r) and C{r) that are now explicitly 
proportional to higher multipole moments of the gravi- 
tational field. According to Eqs. (1.5) and (1.6), then, 
the field at late times will be dominated by a spherical- 
harmonic mode that might be different from the original 
— the final I might be smaller than the initial I, in vio- 
lation of our predictions. 

This, however, is an effect that is clearly nonlinear in 
the gravitational potentials, as the amplitude of the field 



at late times is now found to be proportional to the prod- 
uct of M times Ci or C;, which are already proportional 
to higher multipole moments of these potentials. Our 
calculational methods do not allow us to describe such a 
nonlinear evolution of the scalar field, and we shall have 
to leave this issue unexplored for the time being. But 
while such a nonlinear effect could ultimately be impor- 
tant, we are not aware of any indication in the literature 
that it might already have been seen. Krivan's numeri- 
cal results [41], for example, describe a radiative decay in 
Kerr that is slower, not faster, than what is expected on 
the basis of a linear analysis — he sees an initial Y^ mode 
decaying as t^^ '"', not t~^, while our linear prediction is 

With a modest effort, the work presented in this paper 
could be extended to the more relevant cases of elec- 
tromagnetic and gravitational radiative fields. Methods 
to compute Green's functions for these fields were intro- 
duced in Refs. [43, 44] in the context of a weakly curved 
spacetime with vanishing vector potential A, and these 
methods could easily be generalized to the situation at 
hand. While we can be reasonably sure that these fields 
will give rise to a very similar formulation of the gen- 
eralized falloff theorem, a complete justification of this 
statement must await future work. 



Organization of the paper 

The technical part of the paper begins in Sec. II, 
where we present a detailed description of the station- 
ary, asymptotically-flat, and weakly-curved spacetime in 
which the massless scalar field is propagating; the met- 
ric for this spacetime was already displayed in Eq. (1.2). 
In Sec. Ill we write down the wave equation for the 
scalar field and present its solution in terms of a re- 
tarded Green's function. In Sec. IV we compute this 
Green's function and show that at late times, it reduces 
to the result of Eq. (1.4); some technical aspects of this 
computation are relegated to Appendix A. In Sec. V we 
calculate the late-time behavior of the scalar field and de- 
rive the results of Eqs. (1.5) and (1.6); technical details 
are relegated to Appendix B. In Sec. VI we switch from 
spherical coordinates (r, 9,<j)) to spheroidal coordinates 
(f , 9, (/)) and show that the transformation is responsible 
for the very different description of Eqs. (1.7) and (1.8); 
technical details are relegated to Appendix C. Finally, in 
Sec. VII we give a direct comparison between the two sets 
of initial data obtained in the two coordinate systems. 

Remarks on notation 

Throughout the paper we use geometrized units, in 
which G = c = 1. In the text we often refer to three dif- 
ferent points in spacetime, x = {t,x), x' = {t',x'), and 
x" = {t",x"). We let da = {dt,da) denote partial dif- 
ferentiation with respect to = {t,x'^), da' = {dt>,da') 
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partial differentiation with respect to .t'" = {t',x'°'), and 
da" = {di",da") partial differentiation with respect to 
x"" = (t",^'"^). We let dai3 = d^dp, d^'p' = da'dp', 
da"p" = da"dp", or any other combination indicate re- 
peated differentiation. To generalize this to any number 
of derivatives we introduce a multi- index = 0102 • • ■ an 
that contains a number n of individual (spatial) indices 
— we shall need this only for differentiation with re- 
spect to the spatial coordinates. Thus, On = daia2---a„t 



d 



'N' 



and 9, 



da'' a'' 



II. METRIC OF A WEAKLY CURVED 
SPACETIME 



We consider the weakly curved spacetime of a station- 
ary matter distribution described by a mass density p{x) 
and a mass-current density j{x) = pv, where v is the 
velocity field within the matter distribution. We use the 
bold symbol x to denote the spatial coordinates (x, y, z), 
and a bold symbol like v to denote a vector in ordinary 
three-dimensional flat space; all vectorial operations shall 
refer to this space. 

The metric of such a spacetime can be expressed as 
[45] 

= -{l+2^)dt'^+{l-2^){dx'^+dy^+dz'^)+8iA-dx)dt, 

(2.1) 

where $ and A are gravitational potentials defined by 



^x) 



J \x-x'\ 



(2.2) 



and 



A{x) = - [ d^x'. (2.3) 

J \x — x'\ 

The metric of Eqs. (2.1) (2.3) is a solution to the Einstein 
field equations linearized about flat spacetime; the source 
terms are defined in terms of the matter's stress-energy 
tensor by p = T** and j°- = T*", where a. = 1, 2, 3 repre- 
sents a spatial index. The metric is expressed in the (lin- 
earized) harmonic gauge, and the corresponding gauge 
condition for a stationary situation is V-A = 0. This fol- 
lows automatically from Eq. (2.3) and energy-momentum 
conservation, which produces the (time-independent) 
continuity equation V • j = 0. Throughout the paper 
we will work consistently to first order in the potentials 
$ and A. 

We suppose that the matter distribution is bounded by 
a sphere of radius TZ. In the integrals of Eqs. (2.2) and 
(2.3), therefore, the source point x' is such that r' = \x'\ 
never exceeds TZ. If a gravitational potential is evaluated 
at a field point x such that r = \x\ is greater than TZ, it 
can be expressed as a multipole expansion in powers of 
1/r. Specifically, 



*w=-T-2: 

71=1 



(-1)" 



(2.4) 



and 

n=2 

Here, M = J p{x') d^x' is the total mass of the mat- 
ter distribution, and = J p{x')x'^ d^x' are higher 
multipole moments of the mass density. The multi- 
index N = aia2 ■ - -an includes a number n of spatial 
indices, and we set On = daia2-a„ = da-,da2 ■ ■ ■ da„ 
and x'^ = x'"! x'"^ • • • x'"" . We also have introduced 
J°'^ = J j°-{x')x'^ d^x' as the dipole moment of the cur- 
rent density, and J"'^ = J j°'{x')x'^ d^x' as its higher 
multipole moments. The dipole moment is intimately re- 
lated to J , the total angular momentum of the matter 
distribution; this is defined by J = J x' x j{x') d^x' and 
the relationship is J"^ = —^e'^'"^Jc, with s"-'"' denoting 
the completely antisymmetric permutation symbol. 

The multipole moments and J*^^ are all symmet- 
ric under permutations of individual indices within the 
multi-index N. Because any "pure-trace" quantity of the 
form 5°'^d...a:.b— {^/f) vanishes away from r = 0, only the 
tracefree part of each multipole moment is actually re- 
quired in the expansions of Eqs. (2.4) and (2.5) [46]. For 
our purposes, however, it will not be necessary to explic- 
itly remove the "trace part" of the multipole moments; 
as this does not alter the value of the potentials, we shall 
simply work with the multipole expansions presented in 
Eqs. (2.4) and (2.5). 



III. SCALAR WAVE EQUATION AND 
SOLUTION BY GREEN'S IDENTITY 

We consider a scalar field ^ that satisfies the wave 
equation 



d^{^gg'^f'dp^)=Q, 



(3.1) 



where g is the determinant of the metric. On a spacelike 
hypersurface S we specify initial values for the field and 
its normal derivative. The field's behavior to the future of 
S is then completely determined, and * can be expressed 
as [48] 



G{x,x')dcc''i>{x') 



- ^{x')da'G{x, x') n«' dV; (3.2) 



this is Green's identity in curved spacetime. Here, x de- 
notes a spacetime point to the future of E, and x' is a 
point on the hypersurface, which has n" as its unit nor- 
mal vector and dV as its natural (invariant) volume ele- 
ment. The quantity G{x, x') appearing inside the integral 
is the retarded Green's function for the wave equation of 
Eq. (3.1); this satisfies 



daiV^g'^^dpG) = -4^<54(x - x') 



(3.3) 



7 



together with a condition that G{x, x') vanish if x is in 
the causal past of x'; here 64^{x — x') denotes a four- 
dimensional Dirac distribution. 

For the spacetime of Eq. (2.1), Eq. (3.1) takes the ex- 
plicit form 



(□ + 4$a„ + 8A"ato)* = 0, 



(3.4) 



where □ = —{d/dt)^ + is the flat-spacetime wave 
operator, Eq. (3.3) becomes 

(□ + 4^dtt + 8A''dta)G = -4TTSi{x - x'), (3.5) 

and Eq. (3.2) reduces to 



^{t,x) = G{t,x;0,x')^{0,x')-^{0,x')dt'G{t,x;0,x') [l - 4$(cc')] rfV 

G{t, x; 0, x')da''i'iO, x') - *(0, x')da'G{t, x; 0, a;')] A''{x') (fx', (3.6) 
I 



where ^'(0, x') is the field's initial value on the hypersur- 
face t = 0, and ^{0,x') the initial value for the field's 
time derivative. 

Equation (3.5) can be solved perturbatively about flat 
spacetime. This method of computation originates from 
Ref. [43], and the particular implementation used here 
comes from Ref. [47]; a more detailed exposition can be 
found in Ref. [44]. Wc write 



G{x,x')=G^^\x,x') + G\x,x') + 



where 



G''^\x,x') 



6(1 -t' -\i 



X 



\X — X' 



(3.7) 



(3. 



is the retarded Green's function of the flat-spacetime 
wave equation, □G"''''' — — 47r54(.T — .x'), G^{x,x') is a 
correction of order <I> and A, and the • •" terms repre- 
sent corrections of higher order in the gravitational po- 
tentials. Substituting Eq. (3.8) into Eq. (3.5) reveals that 
the first-order correction to the Green's function satisfies 

UG\x,x') = -4($a« + 2A''dta)G^^\x,x'), (3.9) 

in which the gravitational potentials are evaluated at x. 
The integration of this equation will be carried out in 
Sec. IV. 

The Green's function of Eq. (3.7) can be decom- 
posed into two parts. The first, which is dominated by 
G*''^*(x, x'), is supported entirely on the past null cone of 
the field point x; this is the "direct" part of the Green's 
function [48, 49], which describes the direct propagation 
— at the speed of light — of a disturbance from x' to 
X. The second, which is contained in G^{x,x'), is sup- 
ported inside the past null cone of the field point x; this 
is the "tail" part of the Green's function [48, 49], which 
describes the delayed propagation - at speeds smaller 
than the speed of light — of a disturbance from x to 
x' . This delay is a curved-spacetime phenomenon: it can 
heuristically be understood as arising from a scattering 
of the radiation by the spacetime curvature. Note that 



G^{x,x') contains both a direct part and a tail part; in 
the sequel we shall be interested only in its tail part, 
which we denote G^'''\x,x'). 

Likewise, the scalar field of Eq. (3.6) can be decom- 
posed into direct and tail parts. The direct part of the 
field represents the passage of the main wavefront, that 
is, the arrival at the observation point of that part of 
the initial wave packet that has propagated at the speed 
of light. The tail part of the field is what arrives after 
the main wavefront, delayed by some interaction with 
the spacetime curvature. In the sequel we shall be inter- 
ested only in the tail part of the scalar field, that is, the 
field's behavior after the passage of the main wavefront; 
it is this behavior that reveals the influence of the space- 
time curvature on the radiation. Ignoring terms that are 
quadratic in the gravitational potentials, we find from 
Eq. (3.6) that the tail part of the field is given by 

**^'\t,a;) = /[G*'^"(i,a;;0,a;')*(0,a;') 
A-K J L 

- *(0, x')dt'G^^'\t, x; 0, a;')] (fx', 

(3.10) 

where, as was noted before, G*'"'(a;, x') is the tail part of 
G^{x,x'). 



IV. EVALUATION OF THE GREEN'S 
FUNCTION 

Equation (3.9) has the form of a wave equation with a 
source term, and its solution is [43, 44, 47] 

G\x,x') = G''^\x,x"){^dt'>t" 

+ 2A''dt"a")G^^\x" ,x')d^x" , 

where the gravitational potentials are evaluated at x" . 
Because the flat-spacetime Green's function G^'^^{x" , x'), 
given explicitly by Eq. (3.8), depends only on the differ- 
ence between x" and x' , the derivatives with respect to 
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x" can be converted (by introducing a minus sign) into 
derivatives with respect to x' , and these can be taken 
outside the integral. The first-order correction to the 
Green's function can thus be expressed as 

G^{x,x') = 2dtH'Pix,x')+Adt'a'P''ix,x'), (4.1) 

in terms of two-point functions defined by 

Pix.x') = ^ j G''^\x,x")<^{x")G''^\x\x')(fx" 
2tt J 

= ($) (4.2) 

and 

P"{x,x') = ^ [ G^''*{x,x")A''{x")G^''\x",x')d^x" 

2-K J 

= {A"). (4.3) 

In these integrations, the product of the flat-spacetime 
Green's functions is supported on the two-dimensional 
surface S formed by the intersection of x's past light 
cone and x"s future light cone. The two-point functions 
are therefore the average of the gravitational potentials 
over S, and we have introduced this notation explicitly in 
Eqs. (4.2) and (4.3). We have indicated that the poten- 
tials are functions of the spatial variables x" only; they 
do not depend on the time variable t" . 

It is easy to see that S will have a large area if x 
and x' are widely separated in time. In these circum- 
stances, the spatial vector x" that appears in Eqs. (4.2) 
and (4.3) will be large, and the gravitational potentials 
will be adequately represented by the multipole expan- 
sions of Eqs. (2.4) and (2.5). More precisely stated, in 
Appendix A we prove that if the time delay t — t' satisfies 
the condition 

t - t' > r + 3/ + 2n, (4.4) 

where r = |a;|, r' = |a;'|, and TZ is the maximum radius of 
the matter distribution, then r" = \x"\ > TZ and ^{x") 
and A°'{x") can both be expressed as a multipole expan- 
sion in powers of 1/r". For the remainder of the paper 
we will assume that Eq. (4.4) holds. This means that 
the time t that appears in Eq. (3.10) will be restricted by 
t > r+3£+2TZ, where r is the distance to the observation 
point, C = max(r') the radius of the region that contains 
the field's initial data, and TZ the radius of the region 
occupied by the matter. Notice that we are explicitly as- 
suming that the field's initial configuration is compactly 
supported: the field's initial values are assumed to be 
zero outside of a sphere of radius jC. 

Taking into account Eq. (4.4), we substitute the multi- 
pole expansions of Eqs. (2.4) and (2.5) into Eqs. (4.2) and 
(4.3). Introducing the sequence of two-point functions 

U{x,x') = ^ f G''^\x,x")^G''^\x",x')d^x" 
2n J r 

^ ir-'), (4.5) 
Ua{x,x') = {dar-'), (4.6) 



and generically, 

UN{x,x') = {dNr-^), (4.7) 
we find that Eqs. (4.2) and (4.3) become 

°° (—^)" 

P{x, x') = -MU{x, x')^J2 —T M^Un{x, x') (4.8) 

n=l 

and 

P''(x,x') = r^Ub(x,x') - V r^UN(x,x'). 

(4.9) 

We see from Eqs. (4.1), (4.8), and (4.9) that G'^{x,x') can 
be constructed from the infinite sequence of two-point 
functions Un{x,x'). Our next task is to compute these. 

Suppose that we already have Un{x,x'), and that 
we want U^a, the next member of the sequence. In 
schematic notation we have Un = J GfGdx" and UNa = 
J G{da" f)G dx" , where / = 5jv"(r"-i). By letting 
da" act on the second Green's function we rewrite this 
as UNa = J Gda"ifG)dx" + da' J GfGdx", where we 
have used the fact that the second Green's function de- 
pends only on x" — x'. The first integral is then inte- 
grated by parts, and if we use the fact that the first 
Green's function depends only on a; — x" , we obtain 
UNa ^ daj GfG dx" + da' J GfG dx". The sequence 
therefore obeys the recursion relation 

UNa{x,x') = {da + da')UN{x,x'), (4.10) 

and each member can be obtained from U{x, x') by re- 
peated application of the symmetric derivative operator 

da + da' ■ Explicitly, 

Ua{x,x') = {da + da')U{x,x'), (4.11) 
Un{x,x') = {da,+ da'^)---{da„+da'JU{x, X'),{A.12) 

and we see that the key to obtaining G^{x, x') is to com- 
pute the two-point function U (x, x') defined by Eq. (4.5). 

An almost identical set of manipulations would reveal 
that the two-point functions depend only on t — t' , since 
{dt+dt')UN{x, x') = 0. This property extends to P{x, x') 
and P'^{x,x'), as follows directly from Eqs. (4.2), (4.3) 
and the fact that the gravitational potentials do not de- 
pend on time. Finally, Eq. (4.1) implies that G^{x,x') 
also depends only ont—t'; this conclusion follows directly 
from the stationary property of the spacetime. 

We must now evaluate U{x,x'). In Appendix A we 
show that iov t — t' > r + r' — which is implied by 
Eq. (4.4) — we have [43, 44, 47] 

^(^'^') = ;^lnj^^, (4.13) 

where R= \x — x'\ is the Euclidean distance between x 

and x'. We sec that the spatial dependence of U (x, x') is 
contained entirely in R, which is antisymmetric in x and 
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x' . Thus {da + da')R = 0, and from this it follows that 
Ua{x,x') and all other members of the sequence vanish 
for t - t' > r + r' . Then Eqs. (4.8) and (4.9) imply 
P{x,x') = -MU{x,x'), P''{x,x') = 0, and Eq. (4.1) 
gives 



G\x,x') 



8M{t - 1') 
' [{t - t'Y - E?f 



(4.14) 



This result is valid under the condition stated in 

Eq. (4.4). ¥oi t > f + r + 3r' + 27^, therefore, the full 
scalar Green's function G{x,x') is given by the expres- 
sion of Eq. (4.14), except for corrections of higher order 
in the gravitational potentials. This result is remark- 
able because it shows that the late-time behavior of the 
Green's function is governed entirely by the spherical, 
1/r part of the gravitational potentials. In Appendix A 
we prove that Ua{x,x') and other members of the se- 
quence are nonzero when t — t' < r + r' . This inequality 
is incompatible with Eq. (4.4), which ensures that the 
gravitational potentials can safely be expressed as mul- 
tipole expansions. In other words, Ua{x,x') and other 
members of the sequence are nonzero only where they 
not properly defined. 

V. LATE-TIME BEHAVIOR OF THE FIELD 

In this section we substitute Eq. (4.14) into Eq. (3.10) 
and compute the scalar field '^{t,x). We shall represent 
the field point x in terms of spherical coordinates (r, 9, <j)), 
and the source point x' in terms of (r', 9',(j)'). We shall 
use the short-hand notation = (9, (p) and fi' = {9' , (f)'). 

We need to specify initial conditions for the field and 
its time derivative. We make the specific choices 

*(o, x') = c{T')Yrm, *(o, x') = c{r')Yrm, 

(5.1) 

where Y^{Vl') are spherical-harmonic functions, and 
C(r') and C{r') are two arbitrary functions supported 
in the interval < r' < £ < oo. We therefore sup- 
pose that * is initially a pure multipole field of degree 
I and azimuthal index m (assumed for simplicity to be 
the same for both and ^i). and that the initial data 
is confined to a ball of radius C. Note that C(r') and 
C{r') are two distinct and unrelated functions of r'; the 
overdot simply means that C{r') specifies the initial data 
for the field's time derivative, and it should not suggest 
that C{r') is the time derivative of C{r'). Below we will 
need moments of these functions, Ci and Ci , defined by 



dr' 



(5.2) 



and a very similar relation holding for C;. Notice that 
as defined here, C; is the integral of C(r') multiplied by 
(r')'"""^, not (r')'; the two additional powers of r' are con- 
tributed by the volume element in spherical coordinates: 
dPx' = r'2 dr'dn', where d^' = sin 6*' d9'd6'. 



We consider the behavior of the scalar field at times t > 
r -|- 3£ -|- 27?., as demanded by the inequality of Eq. (4.4). 
Substituting Eq. (4.14) into Eq. (3.10) we obtain 

^>{t,x)=S{t,x) + ^^S{t,x), (5.3) 

where the superscript ">" indicates that Eq. (4.4) is be- 
ing enforced, and where 

S{t,x) = ^ j G\t,x;0,x')C{r')Yr{n') d^x' , (5.4) 

with a very similar relation defining S{t,x). Here the 
Green's function G^{t, x; 0, x') is given by Eq. (4.14), and 
the overdot on S means that this quantity is constructed 
from C in the same way that S is constructed from C; 
it does not mean that S is the time derivative of S. To 
compute S{t, x) and S{t, x) we shall first carry out the 
angular integrations to obtain 

Q{t,x;r') = ^ J G^{t,x;0,x')Y{"{fl')dn', (5.5) 

and then perform the remaining radial integration: 

S{t,x) = J Q{t,x;r')C{r')r''^dr', (5.6) 

with a very similar relation holding for S{t, x). 

The Green's function of Eq. (4.14) can be expressed as 

G' {t, x;0,x') = - -2 T—W^ m2 ' (5-7) 

(i^ — -I- 2rr cos 7 — r'^)^ 



where 



cos 7 = cos 6 cos 9' + sin 6 sin 9' cos{(p — <p') (5.8) 

gives the angle between x and x'. We use the binomial 
theorem to rewrite Eq. (5.7) as 



G' = - 



SMt 



^ ^ p=0 



2rr 



I \ V 



and substitution into Eq. (5.5) gives 



(cos7)P, 



Q{t,x;r') 



p=0 



2rr' 



— — r 



l^Jlpi^)^ (5-9) 



where 



7p(0) = ^ j icos-f)PYi"'{il')dn'. (5.10) 
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These angular integrations are carried out in Appendix of p, and in particular for p < I, the angular functions 



B, where we show that 

2'(; + 2n)!(/ 



7p(0) 



(2/ + 2n+l)!n! 
a p = I + 2n, where n = 0, 1, 2, • 



(5.11) 



for all other values 



vanish. The infinite sum in Eq. (5.9) therefore starts at 
p — I, and it contains even terms only if I is even, or 
odd terms only if I is odd. After some simple algebra, 
substitution of Eq. (5.11) into Eq. (5.9) yields 



0(t x-r') iM( iy+i (^' + ^)" (^^')'* y^(^.f> (/ + l + 2n)!(/ + n)! (2/ + 1)! / 2rr' 
Lj[t,x,r )-iM[ 1) (2Z + l)!!(t2_r2_^/2y+2^i 2^ + 1), (2^ + i + 2n)! _ ^2 _ 



\ 2n 
Zi2 I ' 



n=0 

which we rewrite as 



(2Z + l)!!(t2-r2)' 

.1,. ' ^'+^f. (/ + l + 2n)!(? + n)! (2^ + 1)! / 2rr' 

I +f2_^2_^/2y A. + i\n\ (2Z + l + 2n)!Vi2-r2-r'2y • ^ ' 



In Eq. (5.12), the essential information appears on the 
first line, and the second line (which includes the infinite 
sum) contributes only an overall factor of H-ei+£2^, with 
El denoting a small quantity of order r'"^ /{t^ 



and £2 one of order rr' /{t^ 
of order £, we have that 



r'2) 

■r'2). Because r' is itself 



r 



£1=0 



£2 



i2 - r2 - £2 



£2 = 



^2 _ ^2 _ £2 

(5.13) 

These small corrections to the leading factor of Eq. (5.12) 
will be of no concern to us. 

Substitution of Eq. (5.12) into Eq. (5.6) gives 

s{t,x) = 4M(-i)'+^||±f|a l'l.y^. yri^) 

X [l + 0(£i)+0(£2^)], (5.14) 

and the same relation also gives S{t, x) if Ci is substi- 
tuted in place of Ci . These results can now be substituted 
into Eq. (5.3), and this yields 



rH 



{21 + 1)!! (t2 - r2)'+2 



X ICi 



2{l + 2) 



l + 0(£l) + 0(£2') 



(5.15) 



This expression is valid for t > r + 3£ + 27?., in which 
TZ marks the boundary of the matter distribution, and 
C the extension of the field's initial configuration. Equa- 
tion (5.15) shows that gcncrically, the field behaves as 
rH/{t^ — r2)'+2, but that its behavior goes to /{t^ — 



r2)'+2 instead if \E'(0, a;') = 0, that is, if the initial data 
is time symmetric. 

If t r the general result of Eq. (5.15) takes the sim- 
pler form 

*(t»r,r,!!) = 4M(-l)wS±||^y,"(!!) 



(5.16) 



which shows that in this limit, the field behaves as 
l^2i+z_ This is the behavior of the scalar field near 

future timelike infinity. Equation (5.16) is valid up to a 

fractional correction of order ji^ . 

To obtain the limiting behavior of the field at future 

null infinity, we introduce the retarded time coordinate 

u = t — r and the advanced time coordinate v = t + r. 

We then take the limit of r^^ as u ^ 00, which gives 

rnv = ^,u,n) = 2M(-l)'+ii^-l_y™(0) 
X (c,.-^-^Ci). (5.17) 



At future null infinity, therefore, r^P behaves as l/u'"'"2. 
Equation (5.17) is valid up to a fractional correction of 
order C? jw^. 

Equation (5.15), and its two limiting cases (5.16) and 
(5.17), give us our new generalized formulation of Price's 
falloff theorem for a scalar field ^ starting with the initial 
data of Eq. (5.1) and propagating in the spacetime of 
Eq. (2.1). It states that for f > r 3£ 27^, the field 
behaves as it would in a spherically symmetric spacetime 
of the same mass M\ deviations from spherical symmetry 
do not affect the field at such late times. 
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It has long been known that the radiative decay of 
a scalar field in an asymptotically-flat spacetime goes 
from a behavior near future timelike infinity to 

a behavior at future null infinity. Our deriva- 

tion offers a very clear explanation for this transition: 
The two different behaviors arise naturally by taking 
different limits of the more general description rVP oc 
r'+H/{e - r2)'+2 = (1 - u/v)'+^{l + w/^^)/(2'+V+2). 

VI. RADIATIVE FALLOFF IN SPHEROIDAL 
COORDINATES 

In the preceding sections we found that the late-time 

behavior of a scalar field in a weakly curved spacetime 
is governed entirely by the spherically-symmetric part of 
the scalar potential and that it is oblivious to the 
vector potential A and the nonspherical part of $. In 
particular, we found that the spacetime's total angular 
momentum plays no role in the description of the radia- 
tive decay. 

A very different conclusion seems to emerge from pre- 
vious studies of the late-time dynamics of a scalar field in 
the spacetime of a Kerr black hole [33-41]. These reveal 
a much more complicated picture in which an initial V™ 
mode gets converted into modes of lower degree. As a 
result of mode conversion, '^{t ^ r) behaves as if 
/ = morZ = m+l,as l/i'+™+^ iil>m-\-2 and l — m\s 
even, or as l/i'+'"+2 if ; > m -I- 2 and I — mis odd — see 
Eq. (1.7). Furthermore, r^{v = oo) behaves as 1/u'"'"^ if 
/ = m or I = m+1, or as 1/u' if I > m+2 — see Eq. (1.8). 
According to Hod [37], who derived these results for the 
first time, the difference in the field's falloff must be at- 
tributed to the angular momentum of the black hole, but 
this interpretation is in direct conflict with our previous 
conclusion. 

In this section we show that the results reviewed in 
the preceding paragraph arc largely an artifact of the 
Boyer-Lindquist coordinates used to represent the Kerr 
metric. We show that the mode conversion has nothing 
to do with the spacetime's angular momentum, and that 
we can reproduce all of the preceding results within the 
approach developed in this paper. There is therefore no 
contradiction between the statements made in Sec. V and 
the results derived by Hod [37]. 

It is known that in the limit M ^ (where M denotes 
the mass of the black hole), the Kerr metric reduces to the 
metric of flat spacetime in spheroidal coordinates (r, 9, (p). 
These are related to a Cartesian frame {x, y, z) by the 
relations 

X = psva.6 cos (j), y = psinO sin (f), z = r cos 9, 

(6.1) 

where p = Vr^ -|- e^, with e denoting the ellipticity of 
the coordinate system (the angular momentum of a Kerr 

black hole is equal to Me). The wcak-fleld limit of the 
Kerr metric is then naturally represented as a deforma- 



tion of the Minkowski metric expressed in spheroidal co- 
ordinates. This suggests that a meaningful comparison 
of results must involve a translation of the computations 
carried out in Sec. V into this coordinate system. [No- 
tice that while we use here the notation (r, 9, 4>) for the 
spheroidal coordinates, these are in fact distinct from the 
spherical coordinates used in Sec. V; we have discarded 
the distinguishing carets used in Sec. I to keep the nota- 
tion simple.] 

This is what we shall do in this Section. The fact that 
the coordinates of Eq. (6.1) are not spherical — even 
though they become increasingly spherical as r grows 
much larger than e — will be seen to be crucially im- 
portant. We will show that the difference in the coordi- 
nate systems, which is measured by the parameter e, is 
at the origin of the difference between the two descrip- 
tions of the falloff behavior. While the description is 
very simple in the original spherical coordinates (because 
the late-time Green's function is insensitive to deviations 
from spherical symmetry), it is much more complicated 
in the spheroidal coordinates (because they obscure the 
underlying spherical symmetry of the problem) . The two 
descriptions, of course, are fully compatible. 

We shall retrace the steps taken in Sec. V and trans- 
late these computations into the spheroidal coordinates 
of Eq. (6.1). Wc suppose once more that the field's ini- 
tial configuration is a pure multipole field of degree I and 
azimuthal index m; we describe it, in the spheroidal co- 
ordinates, by Eq. (5.1). The late-time behavior of the 
field is once more given by Eqs. (5.3) and (5.4), in which 
the Green's function can be substituted from Eq. (4.14). 
The volume element is now (Px' ~ (r'^-l-e^ cos^ 9')dr'dVl' , 
and the quantity E? that appears in the Green's func- 
tion is now given by E? = p^ + p'^ — (-? cos^ 9 — 
2 pp' sin 9 sin 9' cos{4> - </>') - 2rr' cos 6' cos 0' - e'^ cos^ 9' . 
The appropriate replacement for Eq. (5.5) is then 

Q{t,x;r') = G\t,x;0,x')Yrm 

x{r'^ + e^cos'^9')dn', (6.2) 

and Eq. (5.6) must be replaced by 

S{t,x) = j Q{t,x;r')C{r')dr', (6.3) 

with a very similar relation holding for S{t,x). The 
scalar field is then given by Eq. (5.3), or ^'(f:,a;) = 
S{t,x) + dtS{t,x). In the following we will go through 
a sketchy evaluation of Q{t,x;r'); this will be suffi- 
cient to show that our method, when combined with the 
spheroidal coordinates of Eq. (6.1), reproduces the Kerr 
results reviewed at the beginning of this section. 

The Green's function of Eq. (4.14) takes the explicit 
form 
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G'^{t,x\Q,x') = - 



SMt 

rp4 



1 + -ipr sin e sin Q' cos{(j) -(/) ) + — ^^^r^ cos 9' + ^ cos^ 6*' 



where 

T^ = t^-p^-p'^ + e^cos^e. 
Making repeated use of the binomial theorem, we re-express the Green's function as 



=2 
2^2 



G\t,x;0,x') 



oo p q 

EEE (n- 



SMt 



(-l)P(p+ 1)! e2"(2rr'cos6')«-"(2ppy-9 



gp-qQq+n 



(6.4) 



(6.5) 



(6.6) 



where S = sin ^ sin ^' cos (0 — </>') and C = cosO'. After substituting this into Eq. (6.2) and rearranging the sums, we 
obtain 



Q{t,x;r') 



CXD OO b 

EE E 



8Mt 

'^A^A^ a\{2c-h)\{h- c)\ 

0=0 6=0 c={b/2} ^ ' ^ ' 



(-l)"+=(a + c+ 1)! e2(''-=)(2rr' co&Qf''-^{1pp'Y 



r2(a+c) 



7„6(r',n) (6.7) 



where {6/2} stands for 6/2 if h is even, or for (6+ l)/2 if 
h is odd, and 

/„6(r', O) = / (r'2 + e2c2)5''C^y;"(0') rffi' (6.8) 
47r 7 

are angular integrations that generalize Eq. (5.10). 

We evaluate the functions /a(,(r',f2) in Appendix C. 
The most important task for our purposes is to determine 
the smallest values of a and 6 such that the functions do 
not vanish; we denote these a* and 6*, respectively. In 
Appendix C we show that a* = m, and that 6* = if 
Z = m, 6* = 1 if / = rn+1, and b* = Z-m-2if / > m + 2. 
In all cases wc find that Ia*b* {r' , fi) ex Y™{il). 

We can now evaluate Q{t,x;r') near future timelike 
infinity (for t ^ r). Assuming that r and r' are both 
much larger than e, and that r ^ r' ^ we use the approx- 
imations p p' r' , and t^. A typical term 
in the infinite sum of Eq. (6.7) is then proportional to 
(rr')°"^+^'=(cos6')2=-''4b/t2(a+c)+3_ r^j^g dominant term 

comes with a = a*, b = b*, and c = {6*/2}, and we find 
that 

Q(t » ^) ^ tL+L+3+e (cos^)^Fr(^^), (6.9) 

where e = if 6* is even, and e = 1 if 6* is odd. Using 
the previously determined values of b* , we obtain the re- 
sults displayed in Table I. Notice that when Z — m is even, 
the angular dependence is given by Y^{Q), the mode of 
lowest degree compatible with the fixed value of m; when 
Z — m is odd, we have instead cos0Y^{^l) oc Y^_^i(Q), 
which is the lowest accessible mode given that mode con- 
version must preserve the parity of the initial angular 
profile. These results are identical to those reviewed at 
the beginning of the section and displayed in Eq. (1.7). 

Let us now examine the behavior of rQ{t,x;r') at 
future null infinity (in the limit v = t + r oo); 



TABLE I: Behavior of Q{t,x; r') near future timclikc infinity. 



If . . . 


Q{t ^ r) is proportional to . . . 


I = m 


{rr')"'Y^{Q)/t^"'+^ 


I + 1 


(rr')'"+^ cos6iyr(f^)/t^'"+^ 


l>m + 2, I — m even 


(rr')'"yr(^^)/t'"^"'+^ 


I > m + 2, I — m odd 


(rr')'"+' cos6iyr(fi)/t'+'"+^ 



we shall express this as a function of w = t — r. 

Here we use the approximations p ^ r, p' ^ r' , 
and ~ uv. After multiplication by r, we find 
that a typical term in the infinite sum of Eq. (6.7) is 
proportional to rt(rr')""''+^'=(cos 6l)2'=-''/„fc/r2(a+c)+4 ^ 

(r')""''+^'=(cos6i)2^-^/„b/(w°+'^+2t>^-'=). In the limit v ^ 
oo, the only surviving terms are those for which c = b. 
The dominant term is then 

Q{v = ^)^\j;^^{coserY::{n), (6.10) 

where we have set a — a* ~ in and b = b*. Using the 
previously determined values of b* , we obtain the results 
displayed in Table II. Notice that the generic behavior, 
for all values of Z > m + 2, is rQ oc It can be 

checked that repeated action of cos 6 on Y^ (fi) gener- 
ates spherical-harmonic functions of higher degree (see 
Ref. [50], Sec. 12.9). If Z - m is even, (cos 6')'-™-2y™ 
can be expressed as a linear superposition of Y^, Y^^27 
1^_^4, . . . , and ij™2- If ^ — TO is odd, on the other hand, 
(cos0)'~™~2y^ can be expressed as a linear superposi- 
tion of yrn+i, i^,™+3, • • • , and yj™2- ^^'^ generic situation 
at future null infinity therefore corresponds to a super- 
position of modes that all come with the same falloff 
behavior. The results derived here are identical to those 
reviewed at the beginning of the section and displayed in 
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TABLE II: Behavior of rQ{t,x; r') at future null infinity. 



If . . . 


rQ{r 


= 3c) is pr(jpi)rti(jiial to . . . 


I = rn 






l = m+l 


ir'r 


+icos6iyr(!:!)/«™+3 


l>m + 2 




-2(cose)'-™-2yr(n)/w' 



Eq. (1.8). 

The description of the radiative fallofi of a scalar field 
presented in this section is much more complicated than 
what was presented in Sec. V. Since the method of calcu- 
lation is identical in both cases, the complexity encoun- 
tered here clearly originates from our use of spheroidal 
coordinates, which obscure the underlying spherical sym- 
metry of the problem. Since the spheroidal coordinates 
are directly related to the Boyer-Lindquist coordinates 
used to describe the Kerr metric, the analysis presented 
here suggests strongly that the apparent complexity of 
radiative falloff in Kerr must be attributed to the coor- 
dinate system, and not to the fact that the spacetime 
possesses angular momentum and higher multipole mo- 
ments. 



VII. DIRECT COMPARISON BETWEEN 
SPHEROIDAL AND SPHERICAL 
DESCRIPTIONS OF INITIAL DATA 

In this section wc attempt a direct comparison between 
the difi'erent pictures of radiative falloff that emerge after 
using spherical coordinates (r, 6, 4>) on the one hand, and 
spheroidal coordinates (f, 6*, ((>) on the other. Notice that 
as in Sec. I, we now use different notations for the two 
coordinate systems. We shall also use the short hands 
Vl = {9, (p) and Cl = {9, cf)). For the purpose of this com- 
parison we will restrict our attention to the field's falloff 
behavior near future timelike infinity. 

The results derived in Sec. VI sec Table I — 
show that when the field's initial data is of the form 
C{f)Yi"^{Cl) when expressed in spheroidal coordinates, 
the evolution produces a dominant late-time mode that 
is either y^(r2) if Z — m is even, or Y^_^i{Cl) if Z — m is 
odd; mode conversion preserves the parity of the original 
mode. In the generic situation (for I > m + 2), this late- 
time mode comes with a falloff behavior that is either 

if I _ ^ Js ^^^^^ Qj. ^-(l+m+2) if /_ TO is odd. 

These results contradict a naive expectation according to 
which a mode should decay as t-(2™+3) ^ ^hile a y^+i 
mode should decay as t-(2™+5) This expectation comes 
from a description in terms of spherical coordinates, ac- 
cording to which a Y{" mode has a falloff behavior given 
by i-(2'+3). i^ is verified for the special cases I = m and 
I = m + 1, but not for the generic situation. 

To shed some light on this contradiction, let us first 
consider initial data of the form ^'(0, a;) = C{f)Y™{Ct) 
with I — m even. We imagine transforming the right-hand 



side of this expression into spherical coordinates (r, 0, (j)). 
This would bring the field's initial data into the general 
form 

*(0,x) =„C(r)y™(0) + „+2C(r)Fr+2(f^) + ---, (7-1) 

in which the radial functions mCir), m+2C{r), and so on, 
are determined by the transformation. We then isolate 
the dominant mode of lowest degree and substitute it 
into Eq. (5.4) to calculate S{t,x), which gives the field's 
late-time behavior. To compute this we first evaluate 
Q{t, x; r') as in Eq. (5.12), which we rewrite as 

«(*'^;^')cx^(^l + :^j (7.2) 

where we ignore all corrections that come with higher 
powers of r. Substituting this into Eq. (5.6) yields 

b(t, oc r {U) 2^ + 2 - A;)! f2m+3+2fe ' 



where 



m+2k 



k=0 



,C{r)r 



dr 



(7.3) 



(7.4) 



are moments of the relevant radial function. 

Under normal circumstances the sum of Eq. (7.3) 
would be dominated by the term k = and we would 
have S oc mCTOr'"F^(f2)/t^'"+^, which goes according 
to the naive expectation. If, however, the first moment 
mC'rn wcrc to vanish, then S would be found to decay 
faster, and the naive expectation would be invalid. As a 
matter of fact, if a number of moments were to vanish, 
so that 



iCm — mGm+2 — ■ ■ ■ — mCl—i — 0, 



(7.5) 



then the sum of Eq. (7.4) would start at fc = 5 (/ — m) — 1 
and we would have 



5(t,a;)cx„G_2r™yr(^)A' 



l+m+l 



(7.6) 



in agreement with the generic results displayed in Table 
I. We see that reconciliation between the results of Sec. V 

and VI relies on the cancellations of Eq. (7.5): the naive 
expectation is false because the correct number of mo- 
ments happen to vanish. 

Before wc investigate whether such a scenario is just 
too far fetched, let us state the conditions that produce 
a reconciliation when / — m is odd. Here Eq. (7.1) must 
be replaced by 



^{0,x) = „+iC(r)y™+i(0) + „+3C(r)yr+3(") 
and Eq. (7.3) becomes 

m+3 , ,o\\ 

Sit,x)^r--^Y-,,in)Y^-^^^ 



(7.7) 



k=0 

m+lCm+l+2k 
f2m+5+2k ' 



fc!(m-|-3- A;)! 



(7.E 
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where m+iCm+i+2fe arc moments of the radial fimc- 
tion m+iC'{r). Normally the sum would be domi- 
nated by the k = term and we would have 5* oc 

„+lC„+lr™+ly™_n(^)/*^"'"^^ in accordance with the 
naive expectation. But if 



m+lCm-\-l — m-t-l^m-t-3 



,+iG_4 = 0, (7.9) 



then the sum of Eq. (7.4) would start at fc = ^ (/ — m — 
1) — 1 and we would have 



S{t,x) (X „,,+iG_2r"^+'ir+i(f^)A' 



j.l+m+2 



(7.10) 



in agreement with the generic results displayed in Table 
I. 

How do Eqs. (7.5) and (7.9) got enforced, given that 
the radial function C{f), in the original spheroidal coor- 
dinates, is arbitrary? While we cannot give a complete 
answer to this question, we can identify the basic mech- 
anism by examining a specific example. Let us consider 
the initial data 



that this result will extend to higher orders. If so, we 
would conclude that qCq vanishes, and we would recover 
the appropriate special case of Eq. (7.5). 

The lesson learned here is that while C{f) is indeed 
arbitrary, the transformation to spherical coordinates 
is sufficiently special that it always produces functions 
mC{r) and m+iC(r) whose first several moments van- 
ish according to Eqs. (7.5) and (7.9). While we cannot 
describe how this happens in complete generality, the 
previous example gives us a vivid illustration of the re- 
markable mechanism at work. 
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^-(0,3;) = C(f)P4(cos6'), 



(7.11) 



in which the angular function is a Lcgendre polynomial 
of order 4. In this expression we substitute f{x, y, z) ob- 
tained from solving — (.t^ + + — e^)f^ — e^z^ = 0, 
which can be inferred from Eq. (6.1). We also substi- 
tute cos^ = z/f{x,y,z) and thus obtain ^'(0,cc) explic- 
itly in terms of Cartesian coordinates. Then we make 
the replacements x — rsm9cos(f>, y = rsinO sincf), and 
z = r cos 6, and finally express ^1/(0, x) in terms of spher- 
ical coordinates. The result is complicated, and we sim- 
plify it by expanding it in powers of and discarding 
terms of order and higher. When we decompose the 
resulting expression in terms of Legendre polynomials, 
we obtain a series of the form 

*(0,a;) = oC(r)Po(cos0) + 2C(r)P2(cos0) + ---, (7.12) 

which is a special case of Eq. (7.1). The relevant part of 
this expression is the spherically-symmetric mode, and 
we find that 

o^w = -1155(7^ 



315 V r"^ 



+ 



3,4 i0;:5 +i0;:6 (7-i3) 



where a prime indicates differentiation with respect to f 
with the result evaluated at r = r. We then see that 



oCo= / oC(r)r2dr = 0(e"), 



(7.14) 



which follows from Eq. (7.13) by straightforward integra- 
tion by parts — the boundary terms are zero because the 
arbitrary function C{f) has compact support. The fact 
that Eq. (7.13) produces a vanishing moment at order 
and is truly remarkable, and we can perhaps believe 



APPENDIX A: TWO-POINT FUNCTIONS 

Our task in this appendix is to derive the expression 
for U(x,x') presented in Eq. (4.13), and to justify the 
inequality of Eq. (4.4). We recall some notation: The 
spacetime point x is decomposed as (t, a;), and similarly, 
x' = {t',x') and x" = (t",a;"); we set r = \x\, r' = \x'\, 
R = \x — x'\, and At = t — t'. The methods employed 
here originate from Rcf. [43], but this particular imple- 
mentation comes from Rcf. [47]; additional details are 
presented in Ref . [44] . 

To compute U{x,x') we substitute Eq. (3.8) into 
Eq. (4.5) and integrate over dt" to obtain 

„ 1 rS(At-\x-x"\-\x'-x"\) 



(Al) 

To eliminate the 5-function we align the z axis with 
the direction of a; — a;' and set x = {x,y,z), x' = 
{x,y, z — R). We then perform the coordinate transfor- 
mation x"{s,a,P) given by 

x" = a; -h ^ — i?2 sin a cos /?, 



y 



= y + T: \/ - R'^ sin asm (3, ( A2) 



1 1 

z = z—-R+-scosa. 

2 2 



In terms of the new coordinates (s,Q!,/3) we have ja; 



[s — Rcosa), \x' — x"\ = 2(5 + -Rcosa), and 
d^x" = ^{s — Rcosa){s + Rcosa)smadadpds. After 
substituting these results into Eq. (Al) and integrating 
over ds, we obtain 

1 f 1 

U(x,x') = — -. — — sinadadd, (A3) 
^ ' ' 477 J |a;"| ^ ' 



15 



where x" now stands for the vector a;"(At, a, /?). 

To evaluate the integral of Eq. (A3) we set 
x" = T] — i] and express the vectors rj and fj 
in terms of their ellipsoidal coordinates {At, a, (3) 
and (s, a, /?), respectively. Wc have 77^ = 
i\/At2 -i?2 sin a cos /3, 77^ iVA<2 - i?2 sinasin/3, 
7?^ = ^Atcosa, ^Vs^ — sin a cos /? = —a;, 

fj^ = ^Vs^ — -R^ sinasin^ = —y, and 77^ = |scosq; = 
—z + ^R. Prom the last set of equations we infer 
s = r + r' . Wc then invoke the addition theorem in 
ellipsoidal coordinates [51] and express l/|a!"| as 

t — TTt — t 

xQr{s>/R)Yr*{a,0)Yr{a,p), (A4) 

where .s< = min(Af, r + r'), s> = max(At,r + r'), and 
P/" and are associated Legendre functions. Substi- 
tuting this into Eq. (A3) we find that integration over 
the spherical harmonics Yi"^{a, (3) gives zero unless / and 
77T. are both zero. The only relevant Legendre functions 
are therefore Po(tt) = 1 and Qo(<i) = 
Gathering the results, we have that 

if At < r + r', and 

if At > r + r'. This last result was previously quoted in 
Eq. (4.13). 

By substituting Eqs. (A5) and (A6) into Eq. (4.11) we 
find that 

Uaix, x') = -2^-L±^ e(r + r'- At), (A7) 

which implies that Ua{x,x') = for At > r + r' . This 
statement carries over to all other members of the se- 
quence, and we conclude that [/(.'/;, x') is the only member 
that does not vanish in this interval. 

The two-point function U{x,x') is the first member of 
a sequence Un{x, x') that arises as a result of expressing 
the gravitational potentials <^{x") and A{x") as multi- 
pole expansions in powers of l/|a;"|. Such expansions are 
meaningful only if 

\x"\> n, (A8) 

where TZ denotes the maximum radius of the matter dis- 
tribution. The two-point functions, therefore, arc defined 
only if the condition of Eq. (A8) can be met. Since x" , 
as defined by Eq. (A2) with s set equal to At, depends 
on the relative positions of x and x' , wc must examine 
under what circumstances Eq. (A8) can be expected to 



hold. As wc shall sec, Eq. (AS) can be replaced by a se- 
quence of stronger inequalities that lead to Eq. (4.4). In 
other words, Eq. (AS) will be enforced, and the two-point 
functions Un{x,x') will be meaningful, so long as 

At>r + 3r' + 211, (A9) 

which is just a restatement of Eq. (4.4). For concreteness 
we shall assume that r > r'. Notice that if Eq. (A9) 
holds, then At is larger than r -|- r', U{x, x' ) is given by 
Eq. (A6), and all other members of the sequence vanish. 

We set x" = a;' -t- ^ and let th e new vector ^ 
have components = |\/At^ — R"^ sin a cos (3, = 
^VAt"^ - i?2 sin Q' sin/?, and = ^{R + At cos a). This 
vector has a length ^ = |^| = 5 (At -|- i?cosa) that satis- 
fies 

l(At-i?)<^<i(At + i?). (AlO) 

On the other hand, 

\^-r'\<\x"\<^ + r', (All) 
and we also have 

r-r'<R<r + r'. (A12) 

These inequalities will now be used to show that Eq. (AS) 
follows from Eq. (A9). 

Equation (A9) can be re-expressed as At — 2 (r' -1-7?.) > 
r -|- r'. But r + r' > Rhy virtue of Eq. (A12), so we have 
At-2{r'+n) > R, which is equivalent to r'+TZ < ^(At- 
R). But it follows from Eq. (AlO) that i(At - R) < ^, 
and we now have r' + TZ < 01 TZ < ^ — r' . This last 
result implies that ^ — r' > 0, and Eq. (All) states that 
^ — r' < |a;"|. So finally we recover TZ < \x"\, and we 
conclude that Eq. (AS) can indeed be deduced from the 
stronger condition displayed in Eq. (A9). 

APPENDIX B: ANGULAR INTEGRATIONS 

In this Appendix we evaluate the integrals 

Ip{n) = ^ j (cos 7)^^;™ (17') dn' (Bl) 

introduced in Eq. (5.10), where cos 7 = cos 6* cos 6*' -I- 
sin 6* sin cos((/) — (f>'). Wc start by recalling the well- 
known property (sec Ref. [50], Exercise 12.4.6) that 
(cos 7)P can be expanded in terms of Legendre polynomi- 
als Pi (cos 7), with the sum over L including even values 
only if p is even, and odd values only if p is odd; in ei- 
ther case the sum is limited hy L < p. Wc combine this 
with the addition theorem for spherical harmonics (see 
Ref. [50], Sec. 12.8), which states that 

^^^(0087) = ^^ E Yif*(0')Fi^(0). (B2) 

M=-L 
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This identity implies that Pl (cos 7)1^;™ (17') integrated 
over dil' is zero unless L = I. It follows that the only 
relevant part of the decomposition of (0037)^ in Legen- 
dre polynomials is 



(cos 7)^ 



2'(2/ + l)p![l(p + 0]! 

{p + i + iy.[i(p-i)]\ 



Pi(cos7) + • 



(B3) 

and this will actually be present in the decomposition if 
p > I and has the same parity as I. The values of p for 
which /p(ri) is nonzero are therefore limited to the set 
{Z, / + 2, / + 4, • • •}. For these, substitution of Eq. (B2) 
into Eq. (B3), and substitution of that into Eq. (Bl), 
allows the integration over dQ' to be carried out. This 
gives 



7p(Q) 



2V[^(p + 0]! 



(B4) 



(p + i + iy.[iip-i)]\ 

This agrees with Eq. (5.11), where we have set p = l + 2n 



APPENDIX C: MORE ANGULAR 
INTEGRATIONS 

In this appendix we give sketchy evaluations of the 
integrals 

Iab{r', O) = / (r'2 + e'^C^)S''C''Yi"'{n') dQ' (CI) 

introduced in Eq. (6.8); we recall the definitions S = 

sin6'sin6'' cos((/) — 0') and C = cos 9'. We are interested 
mostly in determining the smallest values of a and b 
such that lab does not vanish; we denote these a* and 
h* , respectively. We shall derive the general form of 
la'b'i^' but leave all numerical cocfhcients undeter- 
mined. 

The integral lab will be nonzero if either S'^C^Yi^ or 
ga(jb+2Yrn contains a term proportional to Yq when de- 
composed into spherical harmonics. The action of S is 
given schematically by (sec Ref. [50], Sec. 12.9) 

SYr = sin9e^%)Yl!^^' + {)Yl^^'] 

+ sin^?e-^[()r™+i + 01^7+1], (C2) 

where the empty brackets denote undetermined numeri- 
cal coefficients; acting with S therefore lowers and raises 
m by one unit. On the other hand, the action of C is 
given schematically by (see Ref. [50], Sec. 12.9) 

CYr = ( )Yi^, + ( )Yi^„ (C3) 

and it leaves m unchanged. It is clear that to reduce m 

to zero we need a minimum number m of repeated ap- 
plications of the operator S. We therefore have a* = m, 
and it easy to check that S^Y^^ must have the schematic 
form 

m 

s^Yrm = Y^{n)Y,{)ytm+2s{^') 

-Fterms in r^r4,---,r2™; (C4) 



the presence of Y^^{^) comes from the fact that m ap- 
plications of S generate the factor (sin^e**^)™ oc Y^(p,). 
The terms that do not involve F*^ disappear after inte- 
gration over dfi', and we obtain 

m ^ 

la^b = ir(f^) E( ) / + a'C^)C%'.m+2s{^') dn'. 

(C5) 

Repeated application of C now gives 



bvO 



C'Y 



-m+2s 



s'=0 



■m-b+2s+2s' 



,(0'), (C6) 



and substitution into Eq. (C5) gives 

m „ 

la^b = YZ\n)J2{i)e' / Y,"_^_,_,^,,in')dn' 

s=0 •' 

+ E[( )r" + i)e'] [ Yl>_^_,^,^^,A^')dn'}, 

s'=0 



(C7) 



after some reorganization of the sum over s' and elimi- 
nation of integrals that vanish automatically. 

We are seeking the smallest value of h such that Ia*bt as 
now given by Eq. (C7), does not vanish. Let us consider 
a few special cases. 

If / = TO, then the first integral of Eq. (C7) is nonzero 
when b = 2s — 2, and since b cannot be negative, the 
smallest possible value of 6 is 0. For the second integral 
not to vanish we need b =^ 2s + 2s', and again the smallest 
possible value of 6 is 0. We conclude that 6* = when 
I = TO, and that Ia*b* has the schematic form [( )r'^ + 

HI = TO -I- 1, then the first integral of Eq. (C7) is 
nonzero when 6 = 2,s — 1, and the smallest allowed value 
of b is now 1. For the second integral we need b = 1 + 
2s + 2s', which gives rise to the same minimum value 
of 1. We conclude that b* = 1 when I = m + 1. In 
this case, la'b* has the same schematic form as before, 
[( )r'2 + ()e2]rr(f7). 

If Z = TO + 2, then the first integral of Eq. (C7) is 
nonzero when b = 2s, and the smallest allowed value of b 
is again 0. For the second integral we need b = 2-|-2s-|-2s', 
which gives rise to a smallest allowed value of 2. Since 
this value is larger than the first, we can ignore the second 
integral. We conclude that 6* = when Z = to -|- 2, and 
la^b' = ( )e'Y^m. 

If I = TO -I- 3, then the first integral of Eq. (C7) is 
nonzero when b = 2s-|-l, and the smallest allowed value of 
b is now 1. For the second integral we need b = 3+2s+2s', 
which gives rise to a smallest allowed value of 3. Since 
this value is larger than the first, we can once more ignore 
the second integral and conclude that 6* = 1 when I = 
TO + 3. This gives la'b' = ( )e^>r(^^)- 

This pattern continues for all other values of I. For 
l> m — 2we find that b* = l — m — 2, which gives rise to 
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TABLE III: Smallest values of a and b such that the inte- 
gral labir'jfi) is nonvanishing. The Table also displays the 
schematic form of la'b* (r', il), with empty brackets denoting 
undetermined numerical coefficients that depend on I and m. 



Case 


a* 


h* 


Ia'b'{r',n) 


/ = m 


m 





[( )r'2 + ()e2]yr(fi) 


l = m + l 


m 


1 


[( )r'2 + ()e2]yr(n) 


l>m + 2 


m 


l-m-2 


( )e'lT(«) 



Ia*b* = ( )^'^'^m^^)- results arc summarized in Ta- 
ble III, and they were quoted in the paragraph following 
Eq. (6.8). 
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